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Abstract 



A general theory is developed for describing the nonlinear relaxation of 
spin systems from a strongly nonequilibrium initial state, when, in addition, 
the sample is coupled to a resonator. Such processes are characterized by non- 
linear stochastic differential equations. This makes these strongly nonequi- 
librium processes principally different from the spin relaxation close to an 
equilibrium state, which is represented by linear differential equations. The 
consideration is based on a realistic microscopic Hamiltonian including the 
Zeeman terms, dipole interactions, exchange interactions, and a single-site 
anisotropy. The influence of cross correlations between several spin species is 
investigated. The critically important function of coupling between the spin 
system and a resonant electric circuit is emphasized. The role of all main 
relaxation rates is analyzed. The phenomenon of self-organization of transi- 
tion coherence in spin motion, from the quantum chaotic stage of incoherent 
fluctuations, is thoroughly described. Local spin fluctuations are found to 
be the triggering cause for starting the spin relaxation from an incoherent 
nonequilibrium state. The basic regimes of collective coherent spin relaxation 
are studied. 
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I. INTRODUCTION 



The problem of spin relaxation from a state close to equilibrium has a long history and 
is well studied, being related to the description of spin motion in the vicinity of different 
magnetic resonances. This type of spin relaxation is usually characterized by linear differ- 
ential equations, such as Bloch equations. The theory of spin motion close to equilibrium 
has been expounded in numerous literature, among which it would be possible to mention 
several good books [1-7]. 

Essentially nonlinear spin motion arises if the system is prepared in a strongly nonequi- 
librium initial state, e.g. with magnetization opposite to an external magnetic field, and, 
in addition, is coupled to a resonator. Such nonlinear dynamics are commonly treated by 
the Bloch equations supplemented by the Kirchhoff equation for a resonator electric circuit 
[8-11]. However, the phenomenological Bloch equations do not allow for the elucidation 
of different physical processes involved in the behaviour of the system and are not able to 
describe several, probably the most interesting, self-organized regimes of spin motion, as was 
demonstrated in Refs. [12-14]. Some physical models, based on microscopic spin Hamiltoni- 
ans, have also been considered, whose survey can be found in recent reviews [15,16]. But in 
each of these models one standardly studies only some particular substances and considers 
only a part of spin interactions, mainly secular dipole-dipole interactions, and one takes 
into account only some of the known attenuation processes. At the same time, it is evident 
that taking care of only particular model elements can easily lead to wrong physical conclu- 
sions, since real physical materials always include several different characteristics competing 
with each other. The study of nonlinear spin relaxation is of paramount importance not 
solely owing to its theoretical beauty but also because it can be employed in a variety of 
applications, such as the measurement of materials parameters, ultrafast repolarization of 
solid-state targets, creation of sensitive field detectors, usage in quantum computing and 
others, as is discussed in reviews [15,16]. One of the major possible applications is in achiev- 
ing the regime of superradiant operation by spin masers [13,17-19]. Punctuated nonlinear 
dynamics of spin assemblies can also be a new tool for information processing [20]. 

The aim of the present paper is to develop a general theory of nonlinear spin relaxation, 
being based on a realistic microscopic Hamiltonian including, in addition to the Zeeman 
terms, the main spin interactions, and taking account of the different major mechanisms 
of spin attenuation. By considering just some limited models, it is easy to come to false 
conclusions and to predict fictitious physical effects that by no means can exist in real 
materials. It is only by carefully treating different competing mechanisms that one can 
derive reliable physical implications. 

II. BASIC SPIN HAMILTONIAN 

Keeping in mind the applicability of the theory to a wide class of spin systems, we start 
with a rather general Hamiltonian including the major spin interactions the most often met 
in magnetic materials [1-7,21-23]. Let us consider a solid sample containing N vector spins 
Sj enumerated by the index i — 1, 2, . . . , N. The spin operators S« can represent any particles 
of spin S, starting from S — 1/2 to very high spin values. These can be nuclear or electronic 
spins, as in the standard problems of nuclear or electronic spin resonances [1-7,15]. Magnetic 
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molecules, forming molecular magnets, can possess various spins ranging from S — 1/2 up 
to S — 27/2, as is reviewed in Refs. [16,19,24-26]. Bose-Einstein condensates of dilute 
gases (see reviews [27-30]), being placed in optical lattices can form localized clouds with 
an effective spin per site of order 10 2 or 10 3 . Spin dynamics (mainly linear) is an intensively 
developing field of research, called spintronics [31]. 

The Hamiltonian of a spin system can, generally, be separated into two parts, 

tf = E^ + ^£^> (!) 

i i^j 

the first term being related to individual spins, while the second representing spin interac- 
tions. The single-spin Hamiltonian 

Hi = -fioB ■ S t - D{S*) 2 (2) 

consists of the Zeeman energy and the energy of the single-site magnetic anisotropy. Here 
fi = h'-ys, with 75 being the gyromagnetic ratio of a particle with spin S. For electronic 
spins, fi < 0, while for nuclear spins fi can be either positive or negative. The total 
magnetic field 

B = B e z + (B 1 + H)e x (3) 

contains external longitudinal, B , and transverse, Bi, magnetic fields, and also a feedback 
field H of a resonator, if the sample is coupled to a resonant electric circuit. The anisotropy 
parameter D is positive for an easy-axis anisotropy and negative in the case of an easy-plane 
anisotropy. 

The interaction Hamiltonian 

^ = £<^-J^-S, (4) 

a/3 



includes dipole and exchange interactions. The dipolar tensor is 

r 



D?f = % (<W - K <) , (5) 

where a, f3 = x,y, z and 



II Tjj 

Tij = I ^-ij I i ^ij = 5 ^ij = 



This tensor enjoys the properties 



E^ a = o, E^" = o, (6) 

of which the first is exact and the second one is asymptotically exact for a macroscopic 
sample with a large number of spins N ^> 1. A positive exchange integral corresponds to 
ferromagnetic interactions and negative, to antiferromagnetic interactions. 
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It is convenient to represent the Hamiltonians through the ladder spin operators Sf = 
Sf ± iSf. Then the single-spin term (2) writes as 

Hi = -^B Q St - \ ^{B, + H) (St + S-)-D {S*f ■ (7) 
With the notation 

= D%, b^\ (B* -D$- %Bf) , Clj = \ (Bff - iDV?) , (8) 
the interaction Hamiltonian (4) transforms to 

H i3 = a l3 {s\S) - 1 -S+S-)+ b l3 SlS+ + bl j SrST+ 

(9) 



-\-2cijSi Sj -\- 2cjjiS'j Sj Jij (S^ Sj -\- Si S 



The interaction parameters a i3 - = a 3i , bij = bji, and Cy = Cji are symmetric and have the 
property 



a ij ~ E bij ~ E c 'ij ~ ) 



(10) 



j&i) 



j&i) 



following from Eqs. (6). 

The equations of motion for the spin operators are obtained from the Heisenberg equa- 
tions and the commutation relations 

S^ , Sj = 25ijS ri , Si , S^ = ±5ijS^ . 

In order to represent the evolution equations in a compact form, it is convenient to introduce 
the local fields 

^ = ^ E M/ + dtjSr + Cij S+ + J i3 (s* - S> 



and the effective force 
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2cijSj — - dijSj + 2bijSj~ + ^Si — Sj ) 



There is a characteristic frequency, the Zeeman frequency, which we denote as 

u = 



Then as the equations of motion for the spin operators, we obtain 

D 
Ji 



dH ' =-i(uo + tio)S- + fS* + i?- (SrS^ + StSr) , 



dt 

with its Hermitian conjugate, and 

dSf 



= -\ (f + s- + stf) . 



dt 2 

The following description of spin dynamics will be based on these equations. 



(11) 



(12) 



(13) 



(14) 



(15) 
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III. TRIGGERING SPIN FLUCTUATIONS 



Suppose that the spin system is prepared in a strongly nonequilibrium state, being 
polarized along the z-axis. What then could be the triggering mechanisms initiating spin 
motion and their relaxation to an equilibrium state? It is evident that imposing transverse 
magnetic fields would push the spins to move. But assume that there are no transverse 
magnetic fields at the initial time and no transverse coherence is imposed on the system. 
What then would initiate the spin motion? Here it is important to stress the role of local 
spin waves as of the triggering mechanism for starting the spin relaxation. 

The appearance of spin waves is due to the local fields (11). In order to consider spin 
waves, or more generally, spin fluctuations that arise in a state which is not necessarily 
equilibrium, it is appropriate to work with the operator equations (14) and (15). Let us 
define the operator deviation 

5S? = S9- <S?> (16) 

from an average < S" >, which is not necessarily an equilibrium average, but which can 
be an average over a nonequilibrium statistical operator, though such that < Sf > weakly 
depends on the index i, because of which it can be taken out of the sums in Eqs. (11). 
Then, owing to Eqs. (10), we have 



% 



J2 UcijSS] - \ a l3 5Sj + 2bijSSf + Jij (5Sr - 5S~)\ , (17) 



which demonstrates that these local fields really correspond to local spin fluctuations. 

To emphasize the role of the spin fluctuations, let us set B 1 — H — 0, that is, looking 
at the case when the transverse fields do not initiate the spin motion. And, respectively, let 
< Sf >= 0, but the longitudinal polarization be finite, < Sf >^ 0. Then Sf = SSf. The 
behaviour of spin fluctuations is characterized by linearizing Eqs. (14) and (15) with respect 
to the operator deviations (16). The linearization of the single-site anisotropy term in Eq. 
(14) has to be done so that to satisfy the known exact relations for S — 1/2 and S — > oo, 
which can be represented [16,32] as 

SrS? + SIS' = (2 - I) < St > Sr . (18) 



S 

Introduce the single-site anisotropy frequency 



u D = (2S-l)j (19) 



and the effective spin frequency 



< S z > 

uj s = uj - u D ^ — , (20) 
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where u is defined in Eq. (13). Then, linearizing Eqs. (14) and (15), we find 

j t S~ = -iu s S-+ < SI > i , ^8S? = 0. (21) 

The second of these equations, under the initial condition SSf(0) = 0, gives 5Sf = 0. 
Now let us employ the Fourier transforms for the interactions 

a a = T? 2 °fc e * kr<3 ' , a k = <>■:<< ' k r , 

k j&i) 

with the analogous transforms for and Jjj, and for the spin operators 
Using the notation 

a fc = u s + \ ( ^ + J k - Jo) < SI > , (3 k = \ b k < St > , (22) 



ft V 2 ' * V 1 ' ™ ~ h 

from the first of Eqs. (14), we obtain 

j t S k = -ia k S^ + ifaSt . (23) 
Looking for the solution of the latter equation in the form 

S k = u k e-^ + vie™** , 
we find the spectrum of spin waves 

wjb = y/a? k -W . (24) 

In the long- wave limit, one gets 



U! k ~ \U S 



(25) 



where k — > 0, and the summation is over the nearest neighbours. 

In this way, in the spin system there are always transverse fluctuations, which can be 
named spin waves. The latter, as they have been described, are not necessarily the spin 
waves in an equilibrium state, as they are usually understood [33], but are to be considered 
in a generalized sense. Under spin waves, we mean here just transverse spin fluctuations. It is 
these transverse fluctuations that are responsible for triggering the initial motion of polarized 
spins, when there are no external transverse magnetic fields. This is why these transverse 
spin fluctuations can be called triggering spin waves. Taking into account such quantum 
spin fluctuations makes it possible to describe the dynamical regimes of spin motion, which 
do not exist for classical Bloch equations. And it becomes possible to develop a detailed 
picture of how the transverse spin coherence arises from initially chaotic fluctuations. This 
self-organized process of coherence emerging from chaos is one of the most interesting and 
challenging problems of spin dynamics. 
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IV. SPIN EVOLUTION EQUATIONS 



The equations of motion (14) and (15) for spin operators are highly nonlinear. The 
nonlinearity comes from two sources. One is caused by the spin interactions accumulated 
in the local fluctuating fields (11). Another kind of nonlinearity enters through the effective 
force (12) containing feedback fields included in the term H. The treatment of the nonlinear 
spin dynamics will be done here by means of the scale separation approach [11-15,34], which 
is a generalization of the averaging technique [35] to stochastic differential equations. 

Notice, first of all, that there are two different spatial scales. One of them is related to 
local fields (11) describing random spin fluctuations (17), which is characterized by a spatial 
length of the order of the mean interparticle distance ao- At this length scale, chaotic quan- 
tum spin fluctuations prevail. Another length scale is the wavelength A ^> ao corresponding 
to coherent effects associated with the characteristic spin rotation frequency uj s . At the lat- 
ter scale, coherent spin correlations are important. These two different length scales allow 
us to distinguish two types of operators. One type are the local fluctuating fields (11), that 
is, the variables £ > an d £ + , and another type are the spin operators S~ , Sf, and Sf. 
The former, responsible for local short-range fluctuations, can be represented by random 
variables [2,5,11,16,36], while the latter keep track of long-range coherent effects. Respec- 
tively, it is convenient to define two sorts of averaging with respect to the corresponding 
variables. Then the statistical averaging over spin operators will be denoted by the single 
angle brackets < ... >, while the averaging over the random local fields will be denoted by 
the double angle brackets <C ... The latter, treating the chaotic local spin fluctuations 
as white noise, are defined as 

< &>(*) > = < £(t) > = , < UtMt') >= ^ 3 6(t - , 

« &(f)£(0 » = « » = , « C(tW) » = 2 73 5(t - , (26) 

where 73 is the width of inhomogeneous dynamic broadening. 

It is worth stressing that the white-noise approximation (26) is not principal and could 
be generalized to taking into account a coloured noise by including finite relaxation times. 
This, however, would result in much more complicated and cumbersome equations. It is 
therefore more convenient, following the ideas of the scale separation approach [11-15], to 
separate in the temporal behaviour of spin correlations two parts, fast and slow. The fast 
part is connected to the local spin fluctuations described by the spectrum of local spin waves 
(24). The characteristic frequencies of these fluctuations are defined by the near-neighbour 
spin coupling as well as by the applied external magnetic field. Here and in what follows, 
we assume that this external field is sufficiently strong, so that the fluctuation spectrum 
(24) is characterized by the frequencies of the order of the Zeeman frequency ljq, which is 
essentially larger than the frequency terms due to spin interactions. With the time 2tt/ujo 
being the shortest among all other characteristic times, the related fast spin fluctuations can 
be effectively treated as white noise, as is done in Eq. (26). The influence of spin correlations 
slowly decaying in time can be appropriately included into the transverse relaxation time T 2 
determined by the strength of the spin-spin coupling allowing for dipolar as well as exchange 
interactions. This effective relaxation time will also be taken into account in the following 
consideration, together with the effect of line narrowing due to high spin polarization [6] . 
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Averaging over spin operators, because of their long-range role, one can employ the 
decoupling 

<S?S;> = <S?><Sf> (i^j). (27) 

Though this looks like a mean-field approximation, one should not forget that the restricted 
averaging, denoted by the single angle brackets < ... >, by definition, involves only the 
spin degrees of freedom, without touching the stochastic variables £ and £*. Therefore 
the quantum fluctuations are not lost in decoupling (27) but are preserved because of the 
dependence of the spin averages < S" > on the random variables £o and £. Then decoupling 
(27) is termed the stochastic mean-field approximation [11-16]. 

A special care is to be taken in considering the single-site term of Eq. (14). When 
averaging the latter, one has to preserve the exact limiting properties known for S — 1/2 
and S — > oo. The corresponding decoupling, correctly interpolating between the exact 
limiting behaviours [16,19,32] is 

< SrS? + S?Sr > = (2 - i) < S~ >< SI > . (28) 

Thus, for S — 1/2, expression (28) becomes zero, as it should be, and for S — > 00, one has 
2 < S~ >< S? >, again in agreement with the correct asymptotic behaviour. 

Let us average the equations of motion (14) and (15) over the spin degrees of freedom, 
not touching the fluctuating random fields £0 and £. Our aim is to obtain the evolution 
equations for the following variables: The transition function 



1 N 



u 



E < > , (29) 



SN ^ 

describing the average rotation of transverse spin components; the coherence intensity 

1 N 

T - \\ 



showing the level of coherence in the spin motion, and the spin polarization 



1 N 



s = 



SN 



E < s' > . (31) 



defining the average polarization per particle. 

In order to have the evolution equations representing realistic spin systems, but not 
just some unreasonable models, an accurate account must be taken of the main relaxation 
mechanisms. Being based on unrealistic models, omitting important existing attenuation 
processes, it would be easy to fall into the sin of predicting physical effects that in reality 
can never occur. We shall consider the following basic relaxation rates. 

(1) Spin-lattice longitudinal attenuation 71, caused by spin-lattice interactions. The 
corresponding longitudinal relaxation time is Ti = I/71. For different materials, 71 can 
be of different order. At low temperature, when spin-phonon interactions are suppressed, 
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the parameter 71 can be rather small. For instance, in polarized nuclear targets [16] at 
temperature of 1 K, one has 71 ~ s _1 . In molecular crystals below the blocking 

temperature of the order of 1 K, the spin-lattice rate can be between 7! ~ 10~ 7 and 10 -5 
s _1 (see more details in Refs. [16,24-26]). Being small, this relaxation parameter may not 
play an essential role at the initial stage of spin motion, however, it always plays a principal 
role at the late stages of spin relaxation. 

(2) Polarization pump rate 7*, which is added to 71 when the sample is subject to 
a permanent pump supporting a nonequilibrium level of the longitudinal spin polarization. 
This rate can be made much larger than 71. Thus, by means of dynamic nuclear polarization, 
the pump rate for nuclear spins in solids can be as large as 7J 1 ~ 0.01 and 10 s _1 [16]. The 
sum of 7! and 7^ will be denoted as 

r 1 = 7i + 7j. (32) 

(3) Spin dephasing rate 72, due to spin-spin interactions. This rate has been calculated 
by many authors, and the generally accepted value [1-7] writes as 

72 = n p & ^/S(S + 1) , (33) 

where p = N/V is density and n is a coefficient approximately equal to the number of 
nearest neighbours. The process of spin dephasing is mainly due to dipolar forces. Exchange 
interactions slightly narrow the line width (33), yielding [4,6,21] a factor of about 0.8. The 
coefficient in Eq. (33) also depends on the type of lattice, so that the numerical factor here 
is approximate. The value of 72 is usually larger than that of 71. For example, in polarized 
solid targets [16] 72 ~ 10 5 s _1 , in molecular magnets [16,24] it is 72 ~ 10 10 s -1 . Inverse of 
72 defines the spin dephasing time T 2 = 1/72- 

(4) Effective homogeneous broadening 72 (s) takes into account a correction to the spin 
dephasing rate 72, appearing in the case of strongly polarized spin systems. Such a strong 
polarization can be achieved in magnetically ordered materials, by applying strong longi- 
tudinal magnetic fields, or by dynamic polarization techniques. This effective broadening 
reads as 

72 ( s ) = 72 (0)(l- s 2 ) , 72(0)^72, (34) 

where s is an average spin polarization (31) and 72 is given by Eq. (33). The derivation 
of Eq. (34) is explained in Appendix A. Under weak polarization, when s 2 < 1, one has 

72 (s) ~ 72. 

(5) Static inhomogeneous broadening 72 is due to various magnetic defects, crystalline 
defects, field gradients, and a variety of additional interactions always present in any real 
materials [1-7,21,31]. Very often the inhomogeneity develops in matter not because of 
externally incorporated defects, but being due to the internal properties, when a heteroge- 
neous state is more thermodynamically stable than a homogeneous state [37,38]. This, e.g., 
happens in many colossal-magnetoresistance materials [39-41] and in high-temperature su- 
perconductors [42-46], where there appears mesoscopic phase separation. In general, 72 can 
be both smaller as well as larger than 72. However in the majority of cases, to a very good 
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approximation 7^ ~ 72. Summarizing the homogeneous and inhomogeneous mechanisms, 
discussed above, we denote the overall transverse relaxation rate as 

r 2 = 72 (i_ s 2 ) +72 *. (35) 

(6) Dynamic inhomogeneous broadening 73 is caused by fast dynamic spin fluctuations, or 
the local spin waves, discussed in Sec. III. It comes into play through the stochastic averaging 
(26). The value of the broadening, due to local spin waves, is of the order or smaller than 
7 2 [14-16,21]. As is emphasized in Sec. Ill, this dynamic broadening is crucially important 
at the initial stage of spin relaxation, when there are no applied transverse fields. 

(7) Cross relaxation rates arise when there are several spin species in the system. For 
example, if there are two types of spins, S and F, then the dynamic broadening for spin S 
becomes 



73 



yfn/ss + n/sF ■ ( 36 ) 



Cross correlations can influence other relaxation rates, especially if the Zeeman frequencies 
of the spins S and F are close to each other [1-7,15,16]. 

(8) Spin radiation rate ^ r arises when there exist the so-called wave packets of strongly 
correlated spins interacting with each other through the common radiation field. The possi- 
bility of the appearance of such an electromagnetic friction was, first, noticed by Ginzburg 
[47] and later discussed by many authors (see e.g. [48]). This collective radiation rate is 

7r = J- P&S{kL s f , (37) 

where k is the wave vector of the radiating field and L s is an effective linear size of a spin 
packet radiating coherently. Rate (37) has earlier been obtained [47,48] in the classical 
approximation. In Appendix B, we briefly sketch how this rate can be derived in a fully 
quantum-mechanical picture. It is important to stress that the existence of rate (37) presup- 
poses the occurrence of monochromatic radiation with a well-defined constant spin frequency 
u s and wave vector k, and that the radiation wavelength is much larger than the linear size 
L s of a spin packet, so that 

kL s <^l ( k = ^) ■ ( 38 ) 

If these conditions do not hold, no noticeable relaxation rate arises. And under the validity 
of these conditions, one has 

^ w 0.1(A;L S ) 3 < 1 . (39) 

72 

The rate 7r is so much smaller than 72, and usually much smaller than 7I, that it can be safely 
neglected, being absolutely unable to influence the motion of spins. Actually, Bloembergen 
[1] has already analysed this problem and come to the conclusion that the interaction of 
spins through the magnetodipole radiation field is completely negligible. However, one may 
put the following question. Suppose that the considered sample is ideally homogeneous, so 
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that 72 is very small, and let the initial spin polarization be very high, such that Sq ~ 1. 
Then the effective transverse rate (35) at the initial time t — can become rather small. 
Could then the radiation rate (37) play any noticeable role, at least at the very initial stage 
of spin motion? We study this problem below. 

(9) Thermal noise attenuation emerges when the spin system is coupled to a resonant 
electric circuit. The resonator Nyquist noise, due to the thermal fluctuations of current in 
the circuit creates a fluctuational magnetic field, which has to be included in the effective 
force (12). The magnitude of the thermal field, produced by the Nyquist noise, is well known 
[10]. It was found [12-16] that the resulting thermal attenuation is 

7t = ^% coth , (40) 

where r\ is a filling factor, u is the natural frequency of the electric circuit, 7 is the resonator 
ringing width, and ojt = ksT/h is the thermal frequency. Bloembergen and Pound [8] first 
mentioned that, because of the macroscopic number of spins N entering the denominator 
of 7t, the latter is unable to influence any spin motion in a macroscopic sample. This 
conclusion was confirmed by accurate calculations [12-16]. 

(10) Resonator relaxation rate arises when the sample is coupled to a resonant electric 
circuit. Then in the effective force (12) the magnetic field H is the resonator feedback field. 
The role of this field will be thoroughly studied in what follows. 

Summarizing all said above, for the spin averages (29) to (31), we obtain the evolution 
equations 

du 

— = -i(u s + Z -iT 2 )u + fs, (41) 
^ = -2T 2 w + (u*f + f*u) s , (42) 

% = -\Wf + r»)-ri(s-Q, (43) 

supplemented by the initial conditions 

u(0) = u , w(0) = w , s(0) = s • 

In these equations, ( is a stationary spin polarization, the characteristic spin frequency is 

lo s = uj - uj d s , (44) 

with uj given by Eq. (13) and ujo, by Eq. (19). The total longitudinal rate I\ is defined in 
Eq. (32) and the total transverse rate T 2 , in Eq. (35). The effective force is 

f = -^fM>(B 1 + H)+S + lr u, (45) 

where the last term is the friction force due to the interaction through magnetodipole radia- 
tion, and 7 r is the magnetodipole radiation rate (37). Equations (41) to (43) are stochastic 
differential equations, since they contain the random variables £ an d £, whose stochastic 
averages are given in Eqs. (26). The external transverse field B 1 and the resonator feedback 
field H need yet to be specified. 
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V. RESONATOR FEEDBACK FIELD 



The resonator feedback field H is created by the electric current of the coil surrounding 
the spin sample. We assume that the coil axis is along the axis x. The electric circuit is 
characterized by resistance R, inductance L, and capacity C. The spin sample is inserted 
into a coil of n turns, length /, cross-section area A c , and volume V c = AJ. The electric 
current in the circuit is described by the Kirchhoff equation 

L i + * + h §• (46) 

in which Ef is an electromotive force, if any, and the magnetic flux 

$ = — nA c r]m x , (47) 
c 

where 77 ~ V/V c is a filling factor, is formed by the x-component of the magnetization density 

™* = y £ < Sf > ■ (48) 

i 

The electric current, circulating over the coil, creates a magnetic field 
The circuit natural frequency is 

and the circuit damping is 

1= t = 2L = 2Q> (51) 

where r is called the circuit ringing time and Q = uoL/R is the quality factor. Also, let us 
define the reduced electromotive force 

_ cE f 



e f = . (52) 



Then the Kirchhoff equation (46) can be transformed to the equation 

d -§ + llH + [ H(t>) dt' = 7e/ - Anrj ^ (53) 



for the feedback magnetic field created by the coil. 

The feedback equation (53) can be represented in another equivalent form that proved 
to be very convenient for defining the feedback field [12-15]. For this purpose, we involve 
the method of Laplace transforms and introduce the transfer function 
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G(t) = (cosu't - wnw't) e~ 7 ' , (54) 



where 



<J = \Juj 2 — °' 2 



Thus, we transform the feedback-field equation (53) to the integral representation 

H = J* G(t - t') \~ie f {t') - 4nr]m x {t')} dt' , (55) 

in which 

™<x(t) = ^ Pl^oS ^ (u* + u) . (56) 
Let the resonant part of the reduced electromotive force (52) be 

ef(t) = h 2 cos cut . (57) 

And let us introduce the notation 

* = -M- (58) 

As usual, we assume that all attenuation parameters are much smaller than the charac- 
teristic spin frequency uj s . Then Eq. (55) can be solved by an iteration procedure, which in 
first order gives 



h 

Here the coupling function 



= i(au — a*u*) + 2/? cos cot . (59) 



1 — exp{— i{uj — uj s )t — 7t} 1 — exp{— i{uj + u s )t — jt} 



7 + i(u — uj s ) 7 — i{uj + uj s ) 

describes the coupling of spins with the resonator and the function 



(60) 



(3 = | (l - e"*) (61) 

characterizes the action of the resonator electromotive force on spins. In Eq. (60) the 
notation for the natural spin width 

7o = | VP^lS (62) 

is employed. 

The spin-resonator coupling can be characterized by the dimensionless coupling param- 
eter 

72(7^ + A 2 ) 
13 



in which A = uo — \u s \ is the detuning. As is evident from Eq. (60), an efficient spin-resonator 
coupling is possible only when the detuning from the resonance is small, such that 



|A| 



<1 (A = cu- \u s \) . (64) 



When the resonance is sufficiently sharp, so that |A| < 7, then the coupling function (60) 
reduces to 



a 



912 (l - e"*) . (65) 



Thus, the resonator feedback field H is defined by Eq. (59), in which a is given by Eq. (65) 
and 0, by Eq. (61). 

VI. AVERAGED EVOLUTION EQUATIONS 

The resonator field, defined in Eq. (59), has to be substituted in the effective force 
(45) entering the evolution equations (41) to (43). In Eq. (45), we also need to specify 
the external magnetic field B\. In general, the latter may contain a constant part and an 
alternating term. So, let us take this transverse field in the form 

B\ = ho + hi cos uot . (66) 

In what follows, we shall use the notation 

Hoho fiohi , 

,0 = — , * = -y-- (67) 

Equations (41) to (43) are stochastic differential equations, containing the random variables 
£0 and £ describing local spin fluctuations. In order to derive the evolution equations in 
terms of ordinary differential equations, we have to accomplish the averaging over random 
fluctuations. This can be done by following the scale separation approach [11-16], the usage 
of the stochastic averages (26), and by invoking the known techniques of treating stochastic 
variables [49]. 

Keeping in mind that the attenuation parameters are substantially smaller than the 
characteristic spin frequency u s , we notice from Eqs. (41) to (43) that the function u can 
be classified as fast, being compared with the temporal behaviour of the functions w and s. 
The latter play the role of temporal quasi-invariants with respect to u. 

Fist, we substitute into Eqs. (41) to (43) the effective force (45), the resonator field (59), 
and the transverse magnetic field (66). This results in the equations 

du 

— = —i{uj s + io)u — (r 2 — as — j r s)u + fis — asu* , (68) 

(JjL 

^ = -2(r 2 - as - lr s)w + {v?h + f» s - as (u 2 + (u*) 2 ) , (69) 
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dt 

in which 



^ = _( a + lr )w - \ («•/! + f{u) -T x {s-Q + \a{u 2 + (u*) 2 ) , (70) 



fi = -ivo-2i(vi+P)cosut + £. (71) 



Then we solve Eq. (68) for the fast variable u, keeping the quasi-invariants fixed, which 
yields 

u = Moexp | — [ioj s + T 2 — as — -j r s)t — i J £o(*0 dt'} + 

+s J* f^t') exp + T 2 - as - lr s)(t - t') - % J* f (O dt"} dt' . (72) 

Solution (72) must be substituted in Eqs. (69) and (70) for the slow functions w and s. After 
this, the latter equations have to be averaged over time and over the stochastic variables £o 
and £, again keeping the quasi-invariants fixed. To slightly simplify the resulting equations, 
one can take the initial condition for the transition function u in the real form, such that 
Uq = Uq, which is not principal but just makes the equations less cumbersome. 

To present the resulting equations in a compact form, we introduce the effective attenu- 
ation 

rs = 7s + ^ _ ci + n e + a* + p i 1 " 6 )' (73) 



in which 



r = r 2 + 73 -(a + 7r)s. (74) 
And finally, after the described averaging, we obtain the evolution equations 

^ = -2(r 2 - as - lr s)w + 2T 3 s 2 , (75) 



— = -(« + 7> - r 3 s - rx( s - c) . (76) 

These equations are very general. They include various attenuation processes, described 
in Sec. IV, and take into account transverse constant and alternating fields (66), as well 
as the resonator electromotive force (57) entering through function (61). The resonator 
feedback field is responsible for the appearance of the coupling function (65). Notice that 
the radiation relaxation rate 7 r , defined in Eq. (37), enters everywhere together with the 
spin-resonator coupling a. However their values are drastically different. Since 

^ ~ 0.1^- {kL s f « 1 , 
a uj s 

the value of 7 r is so incomparably smaller than a ~ #7 2 , that it is evident, in the presence 
of a resonator, the rate 7 r must be forgotten. 
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Moreover, even when there is no resonator, so that a — (3 — 0, the radiation rate 7 r plays 
no role, since it is much smaller than 72, 7I, and 73. One might think that 7 r could play 
a role in the following unrealistic case. Let us imagine an absolutely ideal lattice with no 
inhomogeneous broadening, that is, let us set 72 = 0, which is certainly a purely imaginary 
situation. Then, according to Eq. (35), one has T 2 = 72(1 — s 2 ). Assume that the spin 
system is completely polarized, with sq = 1. Hence, at the initial time, T 2 = 0. Could then 
the spin motion be started by the term with 7 r ? The answer is evident: As far as the largest 
terms in both Eqs. (75) and (76) are those containing T 3 , the terms with 7 r are always 
negligible, even if T 2 = 0. Even more, functions (30) and (31), by their definition, satisfy 
the inequality 



Therefore, if one sets s = 1, then w = 0, and the term ^y r w simply disappears from the 
equations. Vice versa, if one sets a noticeable wq ~ 1, then s 2 < 1, and r 2 ~ 72 3> j r . In 
this way, the radiation rate 7 r never plays any role in the spin motion, which is in agreement 
with the estimates by Bloembergen [1]. 

Note that the situation in spin systems is principally different from that happening in 
atomic systems. In the latter, both the linewidth 72 = 2|d| 2 /c 3 /3 as well as the collective 
radiation rate 7 r = (2/3)\d\ 2 k 3 N c , where N c is the number of correlated atoms, forming a 
wave packet, are caused by the same physical process, by the interaction of atoms with their 
radiation field. Hence 7 r /72 — N c 1, which results in the coherentization of the dipole 
transitions. This is possible even if kL ^> 1, but the number of atoms in a partial wave 
packet is N c >> 1, since 7 r /72 = N c ^> 1. Contrary to this, in spin systems the linewidth 72, 
given in Eq. (33), is due to direct dipole-dipole interactions, while the radiation rate (37) 
is a result of the spin interactions with their radiation field. This is why in the latter case, 
one always has 7 r 72, and the radiation rate 7 r plays no part in the motion of spins. 

We may also notice that in the effective attenuation (73) the terms due to the presence of 
a constant transverse field are less important than the terms caused by the local spin fluctu- 
ations and by the alternating transverse fields. Therefore, omitting the terms corresponding 
to a permanent transverse magnetic field, we have 



w + s 2 < 1 . 



(77) 




(78) 



Finally, we obtain the evolution equations 



dw 



2(r 2 - as)w + 2r 3 s 2 , 



(79) 



ds 
~dl 



aw - T 3 s - r\(s - C) , 



(80) 



describing the averaged motion of spins. 
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VII. COHERENCE EMERGING FROM CHAOS 



One of the most intriguing questions is how the spin motion could become coherent if 
initially it was not. This is a particular case of the general physical problem of how coherence 
emerges from chaos. 

Being interested in a self-organized process of arising coherence, let us consider the case, 
when there are no external transverse fields pushing spins, that is v\ — f3 — 0. Then Eq. 
(78) yields T 3 = 73. Assume also that there is no pumping, so that 7* = 0, hence r x =71. 
Under these conditions, the initial spin motion, for the time t such that 

71* < 1 , 72* « 1 , 1st « 1 , (81) 
follows from Eqs. (79) and (80) in the form 

w ~ w + 2 7 3 Sq - 72 (l - Sq + K ) w o t , 

s ~ s - [(71 + 7 3 )s - 71C] t , (82) 
where the inhomogeneity coefficient is introduced, 

k = 72 772 • (83) 

If at the initial time no transverse polarization is imposed on the system, and the initial 
coherence function is zero, wo = 0, nevertheless the coherent spin motion starts developing 
according to the law 

w ~ 2 73 sgf K = 0) , (84) 

provided there is an initial longitudinal polarization s 7^ 0. The initiation of the emerging 
coherent motion is caused by local spin fluctuations creating the effective rate 73. Recall 
that in the Bloch equations coherent motion never appears if it is not imposed by the 
initial conditions. Contrary to this, Eqs. (79) and (80) take into account the local spin 
fluctuations triggering the motion of spins. The second of Eqs. (82), keeping in mind that 
usually 71 73 ~ 72, can be simplified to 

s~so(l-73*)- (85) 

At the initial stage of spin motion, their coherence is yet incipient, and the motion is 
mainly governed by quantum chaotic spin fluctuations. The coherentization of the trans- 
verse motion goes through the resonator feedback field and the growing coupling function 
(65). The quantitative change in the spin motion happens when the coupling function (65) 
becomes so large that the term (r 2 — as) in Eq. (79) goes negative, which means that an 
efficient generation of coherence has started in the system. This is analogous to the begin- 
ning of maser generation [15-19]. The moment of time, when the regime of mainly chaotic 
quantum fluctuations transforms into the regime of predominantly coherent spin motion, 
can be called the chaos time. This time t c is defined by the equality as = T 2 , that is by the 
equation 
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as = 72 (l-s 2 )+72 (t = t c ) 



(86) 



From here, the estimate for the chaos time is 



gs 



t c = r\n f±- 2 , (87) 



gs - 1 + s. 



where r is the resonator ringing time defined in Eq. (51). The regime of chaotic spin 
fluctuations lasts till the chaos time (87), after which the coherent stage of spin motion 
comes into play. As is clear from the above equations, the transformation from the chaotic 
to coherent regime goes as a gradual crossover. Notice that the quantity 1 — Sq + « is positive 
since Sq < 1. Then, in order that the chaos time (87) be positive and finite, the inequality 

gs > I - sl + k > (88) 

must hold. For a strong spin- resonator coupling, when gs ^> 1, the chaos time (87) reduces 
to 



t 



C 



S-(l-4 + K ). (89) 



As is seen, there exists a well defined stage of chaotic spin fluctuations, with a finite chaos 
time t c > 0, after which the coherent regime develops, if qsq > 0. The coupling parameter g 
is defined in Eq. (63), from which it follows that one should have w s s > 0. Assuming that 
the initial spin polarization is positive, s > 0, one gets the requirement that u s > 0. The 
latter, by definition (44), is equivalent to the condition u > uj d s. Moreover, the coupling 
function (65) is obtained under the resonance condition (64), which implies that oj s has to be 
close to the resonator natural frequency uo. There are two ways of preserving the resonance 
condition (64). First, one can impose a sufficiently strong external magnetic field Bo, such 
that the frequency u>o, given by Eq. (13), would be much larger than u>d, defined in Eq. 
(19). This becomes trivial for S — 1/2, when w D = 0. If lj 3> uj d , then it is easy to realize 
the resonance condition (64), with uu s ~ u and slightly varying in time detuning A = u — uo s . 

The second way of keeping the resonance condition (64) is by means of the chirping effect 
[16,19]. This requires to vary in time the external magnetic field B so that to maintain the 
equality 



h 

with a fixed detuning. 



+ ( W + lj d s) = A , (90) 



VIII. COHERENT SPIN RELAXATION 

After the chaos time (87), the motion of spins becomes more and more coherent, being 
collectivized by the resonator feedback field, with the coupling function a reaching the value 
(77 2 . At the transient stage, when t > t c but t Ti, we may neglect the term with 71 in Eq. 
(80). Assuming that there is no pumping, that is 7^ = 0, one has Ti = 7i . Let us continue 
studying the case of the self-organized coherent spin motion, when there are no transverse 
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external fields, so that v\ — (3 — 0, hence T 3 = 73. When the coherence is well developed, 
then the main term in Eq. (79) is the first one, while the term with 73 can be neglected. 
Under these conditions, and using expression (35) for the rate T 2 , Eqs. (79) and (80) reduce 
to the form 

^ = -2 72 (l - s 2 + K - gs) w , (91) 



ds 
di 



-gi2w . 



(92) 



The solution of these equations is explained in Appendix C and it yields 

- 7 Manh^U^. 



912 



9 



Here 



t p = I hp 



(93) 



(94) 



is the pulse time showing the duration of the coherent relaxation occurring as a fast pulse. 
The delay time 



r„ 



t = t c + ^\n 

defines the time of the maximal coherence. The pulse width is given by the relation 



(95) 



2 _ }_ 2 



1 + 



\ 



(96) 



in which 

lg = l2{gs c - 1 - «) 

The boundary values w c and s c are 



W c = W + 2 7 3 Sq - 72 (l - Sq + K ) w o t c , s c = s (l - 73^c) , 



(97) 



(98) 



with the chaos time t c given in Eq. (87). Since we are interested in the self-organized 
collective process, when there is no large transverse polarization imposed on the system at 
the initial time, we may set w <^ Sq. Then Eq. (96) simplifies to 



7p = ig + (gi2) 2 w c 



(99) 



The pulse time (94) reads as 
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It is easy to notice that if the spin- resonator coupling is weak, j < 1, then 7 P ~ 7 9 ~ 72 
and t p ~ T 2 . In that case, no self-organized coherence can arise in the system. 

Collective coherent effects appear in the spin motion only if the pulse time r v is smaller 
than the dephasing time T 2 . The inequality t p < T 2 , according to Eq. (100), requires that 

(gs c - 1 - k) 2 + g 2 w c > 1 . (101) 

Three different regimes can satisfy Eq. (101). 

The regime of collective induction happens when 

gs < 1 + k , g 2 w > 1 . (102) 

Then, as is clear from Eq. (97), one has 7 9 < 0, because of which t < t c . This means 
that there is no a noticeable maximum in the coherence function w, since, by definition, the 
delay time (95) should occur after the chaotic stage, so that t > t c . But the latter implies 
that 7 S > 0. 

The triggered coherent relaxation corresponds to 

gs > 1 + k , < g 2 w < 1 . (103) 

And the purely self- organized coherent relaxation takes place when 

gs > 2 + k , w = . (104) 

In this classification, we keep in mind the inequality ^t c <C 1, owing to which w c ~ wo and 
s c w s . The initial coherence is assumed to be weak, so that w « 1. 
For w Sq, the delay time (95) can be represented as 

= % m 4 fa'-~ '-">' , (105) 
2 g 2 w c 



In the case of the purely self-organized coherent relaxation, for sufficiently large coupling 
and initial polarization, such that gs ^> 1, the delay time (105) reduces to 



to = tc + y In 



73*c 



(106) 



where r p = T 2 /gsQ. From these formulas, one sees that if 73 — * 0, then £ - > 00, and no 
coherent relaxation is possible. This emphasizes the crucial role of the local spin fluctuations, 
whose existence results in the relaxation rate 73. 

At the delay time (95), solutions (93) are given by the expressions 

/ n / l + /t\ 2 , . 1 + k , . 

w(t ) =w c +[s c , s(t ) = . (107) 
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And for t >> to, they exponentially decay to the values 



w ~ 4:w(t ) exp(— 27 p t) , 



K 




exp(-2 7p t) . (108) 

At very large times t ~ T 1; the transient equations (91) and (92) are no longer valid. Then 
one has to return to the full equations (79) and (80). With increasing time, the solutions 
tend to the stationary points defined by the zeros of the right-hand sides of these equations. 
Among the relaxation regimes to the stationary solutions, one is especially interesting, going 
through a long series of coherent pulses. This pulsing coherent relaxation takes place under 
a permanent external pumping described by a large pumping rate 7* ^> 71. Then Ti — 7^. 
If also the coupling parameter is sufficiently large, such that g( ^> 1 and 

73 «1, 



9(ll 

then the fixed point of Eqs. (79) and (80) is given by the expressions 

9(12 + 1*2) V 9(ll) 9 V 9(ll) 

corresponding to a stable focus. The relaxation to the stationary solutions (109) realizes 
through a series of sharp coherent pulses, similar to the form of Eqs. (93), with the temporal 
interval between the pulses asymptotically defined by the separation time 

T sep = . (110) 
'2g(it(i2 + l*) 



The number of the separate coherent pulses can be estimated as N sep = \/^{T sep , which 
gives 



N = 

1 " sep 



g((i2 + 12) 
\ Wit ~ ' 



Such a highly nontrivial relaxation regime occurs only under a strong pumping and a suffi- 
ciently strong coupling with a resonator. 



IX. INFLUENCE OF CROSS CORRELATIONS 

When in the sample, in addition to the studied spins, there are spins of other nature, the 
presence of the latter can certainly influence the dynamics of the former. Let us consider 
the case of two types of coexisting spins, S and F. The total Hamiltonian is the sum 

H = Hs + Hf + Hsf (HI) 

of the Hamiltonians for S-spins, F-spins and their interactions. The Hamiltonian H s of 
^-spins is the same as in Eqs. (1) to (4). Let us accept for the Hamiltonian H F of F-spins 
a similar general form 
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H F = Y,H iF + ^H ijF . (112) 

i i^j 

The single-spin terms are 

H iF = -fj, 0F B • F, - D F (F t z ) 2 , (113) 
with the total magnetic field (3). And the interaction terms are given by 

HijF = D ijF F t F j ~ JijF^i ■ Fj , (114) 

a(5 



with the dipolar tensor 

ij 

Assume that the interactions between the >S-and F-spins are represented by the Hamiltonian 



H SF = J2ASi-Fi + J2J2 <Sf if , (115) 

i i^j a/3 

containing the part of the single-site interactions of intensity A and the part of the dipole 
interactions, with the dipolar tensor 



ij 



In particular, these could be hyperfine interactions between nuclear and electron spins 
[15,50]. 

We employ notation (8) for the interaction parameters of S-spins and an equivalent 
notation for the interaction parameters aij F , hjF, and Cij F of F-spins. Similarly, we define 
the interaction parameters 

= h = \ ~ 4f ~ ^t) , = \ - ( 116 ) 

for the spin cross interactions. 

The local fields (11), acting on S-spins, are generalized to the form 

& = I E M/ + 4 SJ + c tJ S+ + J^St - S£) + a^F/ + c^F" + c-F/ 



i 



£ — t ■ ^ 2>CijSj dijSj + 2b ij Sj' + Jij{S ri Sj )+ 

+2%-F/ - 1 a l3 F~ + •%/•;/ . 
Analogous local fields act on F-spins, 

^.of = tJ2 [ovfF; + c* jF F~ + CljF F+ + J ijF {F* - F/) + a l3 S* + c^Sj + 



(117) 
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i \ 1 

£,f = ^ X! 2°ijFFj — ^ o-ijFFj + 2b ijF F+ + Jij F {F i — Fj ) + 



Instead of one effective force (12), we have now two forces 



f=-l. liQ (B 1 + H) + j:A i Fr + Z, 

fF = -\ ^of(B 1 + H) + % - ASr + £ F . 
In addition to frequency (13), let us introduce the effective frequencies 



(118) 



_ ^OF D 
Uqf = ^— -DO ) 



(119) 



(120) 



The Heisenberg equations of motion for the system with Hamiltonian (111) yield the 
equations for S-spins 



dt 



-i (coo + eF? + Co) S- + S Z J +'-D (SrSf + SfSr) , 
dSf -l(f + Sr + S+f), 



dt 



(121) 



and the equations for F-spins 
dF,r 



dt 



u 0F + eSf + Zof) Ff + Fff F + — Dp (f~F* + F?Fr) , 



dFf 
dt 



= ~\ {ftF- + F?f F ) . 



(122) 



Again we assume that the sample is inserted into the coil of a resonant electric circuit. The 
feedback field acting on the sample is given by Eq. (53) or (55), where now the magnetic- 
moment density is 



TO, 



i=l j=l 



(123) 



with N F being the number of F-spins. 

Averaging Eqs. (121) and (122), we derive the evolution equations for functions (29), 
(30), and (31), corresponding to S-spins, as well as the equations for the functions 



a, = T r T ~ X " F i ' • 



FN 



F i= i 



(124) 
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WF ~ F 2 N F (N F - 1) 



N F 

xTv^i) g < w > • < 125 > 



SF ^g<^>, (126) 

describing F-spins. In this notation, the transverse magnetic-moment density (123) is 

1 1 

"i* = g PPo S ( u * + u ) + 2 PfPofF(u* f + Mp) , 

where is the density of F-spins. 

The analysis of the evolution equations for the combined system of S- and F-spins is the 
same as has been given above for one type of spins S, with the difference that all expressions 
become much more cumbersome. Again it is possible to show that in the triggering of spin 
motion an important role is played by the coupled S — F spin fluctuations, which yield the 
dynamic relaxation rates 73 and 73^ defined by the relations 

li = l 2 ss + l 2 SF, iIf = 1ff + 1fs , ( 127 ) 

where 

Iff *PF^f yjF{F + 1) , lFS ^VPP~F^S. 
The effective frequencies of S- and F-spins, respectively, are 

ujg = ujq — uj£)S + espS , uop — ^of — ^df^f + ^sF , (128) 
where up, is given by Eq. (19) and 

ludf = (2F - 1) ^ . (129) 

We shall not overload this paper by a detailed exposition of various cross correlations 
resulting from the complicated system of the coupled evolution equations for S- and F-spins. 
Let us only emphasize the existence of a rather nontrivial nonlinear effect of mutual spin 
interactions through the resonator feedback field. Calculating the latter from the integral 
representation (55), with the transverse magnetic density (123), and substituting this into 
the evolution equations results in an effective mutual influence of spins through the feedback 
field. If the resonator is tuned to the characteristic frequency ujs of S-spins, then for the 
latter, we derive the evolution equations similar to Eqs. (79) and (80), but with the effective 
spin-resonator coupling 

77cMs (, . PfVof£s f F\ 

9= , 2 , A2 x 1 + > (130) 

72 (T + A 2 ) \ ppoOJp J 

instead of Eq. (63), and with 73 given by Eq. (127). Depending on the spin characteristics, 
coupling (130) can substantially surpass the value of Eq. (63). This is because the subsystem 
of F-spins, coupled to a resonator, becomes itself a kind of an additional resonator for S- 
spins. 
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X. CONCLUSION 



A general theory is developed for describing nonlinear spin relaxation, which occurs 
when the spin system is prepared in a strongly nonequilibrium state and when the sample is 
coupled to a resonator electric circuit. A strongly nonequilibrium initial state can be realized 
by placing a polarized sample into an external magnetic field, whose direction is opposite 
to the sample magnetization. Nonlinearity in spin relaxation comes from direct spin-spin 
interactions and from their effective interactions through the resonator feedback field. Direct 
spin interactions are responsible for the appearance of local spin fluctuations, playing a 
crucial role at the starting stage of relaxation. The resonator feedback field collectivizes the 
spin motion, leading to coherent collective relaxation. The developed theory is based on a 
realistic Hamiltonian containing the main spin interactions. The role of various relaxation 
rates is thoroughly analysed. 

The aim of the present paper has been to develop a general theory providing an accurate 
and realistic description of nonlinear spin relaxation. This theory can be employed for 
a large class of polarized spin materials. Applications to particular substances require a 
special consideration and separate publications. There exists a large variety of materials 
that can be treated by the developed theory. Just to give an example, we may mention the 
class of molecular magnets [16,19,24-26]. For instance, the molecular crystal Vi 5 is made of 
molecules of spin 1/2, so has no magnetic anisotropy. Its nonlinear spin relaxation can be 
realized in a rather weak external field B > 1 G. The molecules Mn 12 and Fe§ possess the 
spin S = 10. They form crystals with density p ~ 10 21 cm' 3 . The anisotropy frequency is 
ujd ~ 10 12 s _1 . At low temperatures below about 1 K, the molecules can be well polarized, 
with the spin-lattice relaxation parameters 71 ~ 10" 5 — 10~ 7 s _1 . The line width is caused 
by rather strong dipole interactions, with 72 ~ 10 10 s _1 . The condition u > uj d can be 
reached for B > 10 5 G. In the molecular magnet, formed by the molecules Mn 6 , whose 
spin is S — 12, the magnetic anisotropy is much weaker, with ujd ~ 10 10 s _1 , being of the 
same order as 72 ~ 10 10 s _1 . Therefore the required magnetic field is not high, Bq > 10 3 
G. Coupling a molecular crystal to a resonant circuit with the natural width 7 = u/2Q, 
where Q is the resonator quality factor, one can attain the values of the coupling parameter 
as large as g ~ Q ~ 10 4 . With such a strong coupling, the influence of the resonator 
feedback field outperforms other relaxation mechanisms, producing fast coherent relaxation, 
with relaxation times r p ~ 10~ 13 s. Such a fast reorientation of the magnetic moment can 
result in the emission of radiation pulses of high intensity. 
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Appendix A: Effective Homogeneous Broadening 



The homogeneous broadening, existing in spin systems, arises from spin-spin interactions 
and is usually expressed through the moments M n , which may depend on the level of the 
longitudinal polarization s, provided the latter is sufficiently large. The moments have been 
calculated in a number of works [1-7,21]. The most general and exact formula, relating the 
effective broadening with the moments, can be found in Abragam and Goldman [6], which 
for the Gaussian line shape is 



ttMI(s) 
2[M 4 (s) - Ws)\ ' 

The Lorentzian line shape yields to practically the same expression, with a slightly different 
coefficient. The broadening 72(5) for the Lorentzian line is -Jit of the Gaussian broadening. 
The dependence of the moments on the polarization has been accurately calculated [6], 
yielding 

M 2 (s) = M 2 (0)(1 - s 2 ) , M 4 (s) = 2.18Mf(0)(l - s 2 )(l - 0.42s 2 ) . 
Substituting this into 72 (s), and taking into account that s 2 < 1, we obtain Eq. (34). 
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Appendix B: Spin Radiation Rate 

To get a fully quantum-mechanical microscopic picture of spin interactions with electro- 
magnetic field they radiate, one has to add to the spin Hamiltonian (1) the field Hamiltonian 



8tt 

where E = E(r, t) is electric field and H = H(r, t) is magnetic field, and the operator energy 
of spin-field interactions 

N 

i=i 

where Hj = H(rj,t). From the Heisenberg equations of motion for the field variables, one 
finds the vector potential 

1 f . ( , |r — r'l \ dr 



A(r,t) = - /j(r',t- 
in which the current density is 



i /i ' 
r — r 



N 

j = -c/io^Si x W(r-ri) • 

The vector potential A, = A(r i: t) can be represented as 

Aj = Aj + Af + Aj , 



where 



with the notation 



/2 = — (e^ - ie^) , /2 = /^e 2 



From here, we get the magnetic field Hj = H(r i? t) acting on an i-th spin as Hj = Vj x Aj, 
which gives the field 

Hj = H~ + H+ + Hj , 

in which 

" /2* - (g ■ njjjnjj g - 3(/2* • n^jnjj / 9 



HT = -E 



ST I * - — 

c 
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If the spins on different sites move independently of each other, so that the single-spin 
terms in the above sums chaotically oscillate, then the average magnetic field acting on 
each spin from the radiation of other spins is zero. Noticeable action of other spins can 
arise only if there exist the groups of spins, the so-called spin packets, which are strongly 
correlated, moving together. A substantial mutual interaction between spins, caused by 
their electromagnetic radiation, can appear only when this radiation is monochromatic, 
with a well-defined spin frequency u s , the related wavelength A = 2nc/uo S) and wave vector 
k = LOgjc. This radiation can collectivize spins in a spin packet of size L s , provided that 

kL s < 1 . 

When the radiation wavelength A is much larger than the system length L, then L s = L. 
This, however, is not compulsory, and the size of a spin packet can be much shorter than 
L, but it should be much larger than the mean interspin distance. Thus, inequality (38) is 
a necessary condition for the appearance of collective effects. 

Under condition (38), the above magnetic fields can be simplified, averaging them over 
spherical angles. The resulting expressions have to be added to the magnetic field in the 
effective force (12), which acquires one more term, being the friction force 

/' = {lr ~ i$u)u , 

in which the collective radiation rate and frequency shift are 

£^ sin( krij) ^ , 
7r = 7oE I e(ct-r tJ ) , 

j KT 'ij 

^ 7o g^(^) e(ct _ ry)> 

j KTij 

where 

is the single-spin natural width, O(-) is a unit-step function, and N s 
of spins in a spin packet. These formulas can be further simplified to 

2 

7r = loN s = — fj, 2 Sk 3 N s 

and 

The frequency shift is very small, even as compared to 72, since 

— -0.1(A;L S ) 2 «1. 

72 

Of course, such a small shift can be omitted, being negligible as compared to 72 and the 
more so as compared to u s . And for the radiation rate 7 r , substituting there N s = pL 3 , we 
obtain Eq. (37). 



= pL z s is the number 
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Appendix C: Transient Stage of Relaxation 

After the chaotic stage of spin fluctuations, the transient stage comes into play, charac- 
terized by Eqs. (91) and (92). The latter, by introducing the function 

y = 72 (l - s 2 + k - gs^j 
and keeping in mind a sufficiently large coupling parameter g ^> s, rearrange to 

dw dy o 

~dt = ~ 2yW ' dt = {9l2) W ■ 

Differentiating the second of these equations, we have 

d 2 y dy 

-d¥ + 2y Tt =0 > 

which yields 

dy |2 2 

with 7 p being an integration parameter. This Riccati equation possesses the solution 



y = 7 P tanh 



in which 7 p r p = 1 and t is another integration constant. Inverting the dependence of y on 
s for s 2 < 1, we get 

y 1 + k 
s = 1 . 

gi2 g 

This gives the second of Eqs. (93), while the first of solutions (93) follows from Eq. (92). The 
integration constants 7 P and to are defined by the initial conditions, which for the transient 
stage are w c = w(t c ) and s c = s(t c ). 
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